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Weyl semimetals are a class of topological materials that exhibit a bulk Hall effect due to time- 
reversal symmetry breaking. We show that for the idealized semi-infinite case, the Casimir force 
between two identical Weyl semimetals is repulsive at short range and attractive at long range. 
Considering plates of finite thickness, we can reduce the size of the long-range attraction even making 
it repulsive for all distances when thin enough. In the thin-film limit, we study the appearance of an 
attractive Casimir force at shorter distances due to the longitudinal conductivity. Magnetic field, 
thickness, and chemical potential provide tunable nobs for this effect, controlling the Casimir force: 
whether it is attractive or repulsive, the magnitude of the effect, and the positions and existence of 
a trap and antitrap. 


PACS numbers: 04.50.-h 11.15.Yc 73.43.-f 78.20.Jq 

In 1948, Casimir [T] showed that quantum fluctua¬ 
tions in the electromagnetic field cause a force between 
two perfectly conducting, electrically neutral objects. 
This has since been extended to other materials 
Throughout this time, Casimir repulsion between two 
materials in vacuum has been a long sought after phe¬ 
nomenon There are principally four categories in 

which repulsion can be achieved: (i) modifying the di¬ 
electric of the intervening medium 011 [ 7 ], (ii) pairing a 
dielectric object and a permeable object 0 (such as with 
metamaterials i):( iii) using different geometries [9Hn], 
and (iv) breaking time-reversal symmetry min. In this 
paper, we are concerned with Casimir repulsion in iden¬ 
tical time-reversal broken systems. Specifically, we will 
study how Weyl semimetals with time-reversal symmetry 
breaking can exhibit Casimir repulsion. The key ingre¬ 
dient to Casimir repulsion in this paper is the existence 
of a nonzero bulk Hall conductance axy ^ 0, a^y = —cTyx 

[Sl¬ 
it is a general theorem that mirror symmetric objects 
without time-reversal symmetry breaking can only at¬ 
tract one another with the Casimir effect m- This is 
understood with the Lifshitz formula 0 where if we have 
two materials characterized by the two reflection matri¬ 
ces Ri and i ?2 and separated by a distance a in a parallel 
plate geometry, we have 

Ec{a)=h j 2A J ^ tiTogp-(1) 

where the trace is a matrix trace and Qz = . 

This integral generally yields an attractive force; how¬ 
ever, if we break time reversal symmetry, obtaining an¬ 
tisymmetric off-diagonal terms in the reflection matrix 
Rxy = —Ryx there is the possibility of Casimir repul¬ 
sion [16] . One candidate is a two-dimensional Hall mate¬ 
rial and similarly, another is a topological insulator 
where the surface states have been gapped by a magnetic 
field [TS] [17] . A Hall conductance does not guarantee re¬ 
pulsion; longitudinal conductance can overwhelm any re¬ 
pulsion from the Hall effect (although the magnetic field 



FIG. 1. The setup we will consider here is two Weyl semimet¬ 
als separated by a distance a in vacuum and with distance 
between Weyl cones 2b in k space (split in the 2 direction). 

can lead to interesting transitions [18]), and a Hall effect 
that is too strong can suppress Kerr rotation and hence 
lead to attraction. The latter case is an interesting phe¬ 
nomenon where “more” of a repulsive material can lead 
to attraction. 

The material we are interested in is marginal in both 
the case of longitudinal conductance and in an over¬ 
whelming Hall effect: Weyl semimetals [14] with the 
Casimir setup seen in Fig. [^ and the resulting normal¬ 
ized Casimir pressure seen in Figs. [2]|^ These materials 
have linearly dispersive band structures characterized by 
Weyl nodes with different chiralities and characterized 
by a chiral anomaly nn. Clean Weyl semimetals at zero 
temperature have a zero dc longitudinal conductivity and 
optical conductivity Re[cra;a:] oc w [5D]. Additionally, they 
exhibit a bulk Hall effect exemplified in the dc limit by an 
axionic field theory m where in addition to the Maxwell 
action, we have 

2 p 

Sa = J (2) 

where 0{r,t) = 2b • r — 26ot and 2b is the distance be¬ 
tween Weyl nodes in k space whereas 2bo is their energy 
offset, e is the charge of an electron, h is Planck’s con¬ 
stant, c is the speed of light, is the electromagnetic 
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field strength tensor, and is the fully antisymmetric 
four-tensor. Inversion symmetry breaking Weyl semimet¬ 
als, on the other hand, do not exhibit a dc Hall effect 
and therefore will not see the effects described in this pa¬ 
per. The electrodynamics of this were investigated in 
Ref. [23] where the authors even comment on the possi¬ 
bility for a repulsive Casimir effect. 

The marginal nature of Weyl semimetals makes them 
prime candidates for tuning the Casimir force between 
attractive and repulsive regimes. In constructed Weyl 
semimetals made of heterostructures of normal and topo¬ 
logical insulators |53| an external magnetic field can con¬ 
trol the Hall effect and hence the repulsive effects. 
Additionally, some of the first materials that have been 
predicted were pyrochlore iridates [25]; these could also 
see a repulsion tunable with carrier doping or an addi¬ 
tional magnetic field. 

In a real material and experiment at finite tempera¬ 
tures, disorder and interactions should be taken into ac¬ 
count, and in Weyl semimetals they lead to a finite dc 
conductivity OEllEil. We simulate this effect in the 
latter part of this paper by raising the chemical poten¬ 
tial in our clean system, leading to intraband transitions 
that contribute to the longitudinal conductivity (in the 
dc limit these are singular contributions). 

To begin, we consider two semi-inhnite slabs of a Weyl 
semimetal (z < 0 and z > a to be precise), neglecting all 
frequency dependence to the conductivities by assuming 
the electromagnetic response is captured by Eq. ([^. The 
result is just a material that is solely a bulk Hall material 
with current responses given by the Hall conductivity 
o'xy = e^&/27r^/i. This response can be encoded in the 
dielectric function so that Exx = ^yy = ^zz = 1 and Cxy = 
~^yx = io'xy!^- With this set up, if an incident wave k 
hits such a material it will break up into two different 
polarizations in the material that satisfies = kx, 
= ky, and (kf)'^ = kz{kz ± axylc). Additionally, 
the two elliptical polarizations Dj- = e(a;)Ej- are Dj- oc 
-^{kz ± o'xylc)ei =F ikfe 2 where 62 is perpendicular to 
the plane of incidence and ei = 02 x k^. Notice that for 
kz < (^xyjc^ one of the polarizations is evanescent. 

The incident and reflected polarizations can be bro¬ 
ken up into transverse electric (TE) and transverse mag¬ 
netic (TM) modes, and the reflection matrix i?(a;,k) 
just connects incident to reflected = 

i?(a;, k)(if(]’'^, A™)^. As shown in the Lifshitz formula 
Eq. 0, we need the imaginary frequency reflection ma¬ 
trix. If we let uj ^ iuj and dehne + k"^ + ky, 

the reflection matrix for a semi-infinite slab of this bulk 
Hall material is 


Roo i^Qz) 


^ f Q— ^xyj^ Q+ T ^Q_z \ 
(^xyjc \ Q+ ~ 2^2 Q— ~ ETxylcj 


( 3 ) 


where we have defined for brevity Q± = 
^2qz{^q1 + (^xylc^ ± 92 ) (the real frequency ver¬ 
sion of i?oo is found in the Supplemental Material [28] 1. 
Inspecting Roo{i<lz), we see that the reflection matrix 



FIG. 2. (Color online) The normalized Casimir force between 
two semi-infinite bulk Hall materials. Repulsion is seen for 
Oxyajc < 4.00. Po is the distant-dependent ideal Casimir 
force [Tj. For Oxyajc —>■ 00 , PdPo —^ 1. 


only depends on the ratio cqzluxy This dependence 
has implications for the Casimir force. After changing 
variables to solely we can inspect the Casimir 
pressure Pc{a) = —dEc/da, and we have an expression 


Pa = 


2hc 


I dqz ql 


1,/c 


,2qza 


with a function 


g(^^^Pj^,2qza) written out in the Supplemental Material 
P5] for completeness. If we then change variables to 
X = 2aqz and normalize by Casimir’s original result for 

t 2 

perfect conductors Pq = — [I|, we can write the 

equation for the pressure as Pc/Po = ficTxya/c). 

With this formulation, we plot normalized force Pc/P q 
as a function of axyd/c obtaining the single function seen 
in Fig. We see that for Uxya/c < 4.00 we obtain 
repulsion whereas for OxyO^/c > 4.00 we obtain attrac¬ 
tion. Thus, these similar materials trap each other at a 
fixed distance simply dependent on the Hall conductivity, 
OTrap ~ mscrt the value of Oxy = e^b/2'K'^fi 

into this expression, we find oxrap ~ 860/6. This means 
that if 1/6 ^ 0 (nm), then oxrap ~ 0 (^m) quite reason¬ 
able. 


As the distance between the materials gets long, 
Pc/Po 1. This behavior is markedly different from 
the thin film Hall case obtained by Tse and MacDonald 
in Ref. [Hj and Rodriguez-Lopez and Grushin Ref. HZ). 
They found a small [two-dimensional (2D)] quantum Hall 
effect implies a quantized and repulsive Casimir force at 
long distances. In our case, we get attraction at long 
distances for a bulk Hall material independent of the 
magnitude of the Hall effect. To resolve this seeming 
inconsistency, imagine a finite thickness of the bulk Hall 
material of thickness d, then the two-dimensional con¬ 
ductivity axyd diverges as d ^ 00 , and in the case of a 
2D quantum Hall plate with inhnite Hall conductivity, 
the Casimir effect is attractive and approaches Pq. 

To make this argument more precise, one can actually 
find the reflection matrix for a bulk Hall system of thick¬ 
ness d and the result is (derivation of Rd depends only on 
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the axionic action Eq. ([^ and can be found in the Sup¬ 


plemental Material [55] j calculated for real frequencies) 


Rd{iqz) 


( Rxx 

\ — Rxy 



(4) 


with 


Rxx = - sinhQ+d -|- ^ coshQ+d - Q+ sinQ.d - cosQ-d)/D, (5) 

Rxy = — 5 ^(Q+ sinhQ+d -|- 2qz coshQ+d — Q- sinQ-d — 2qz cosQ-d)/D, ( 6 ) 


where 

D = (Q+ -I- 5 ^) coshQ+d -|- {2qzQ+ + ^Q-) sinhQ+d -I- (Q?. - 5 ^) cosQ-d -|- {2qzQ- - ^Q+) sinQ^d. (7) 


It can be shown that limd_>.oo Rd(iqz) = Rao(iqz)- 

Similarly, in the limit of d —> 0, if we keep = cr^yd 
constant, we obtain what was found in Ref. 


lim Rd{iqz) 

d —^0 


1 

l + {al^/2cy{ ofj2c 


/ 2 c \ 

-{afj2cf) • 


For the rest of our discussion, define Roiiqz) = 
Rd{iqz) with = cr^yd held constant. 

With the correct limits identified, we first notice that 
we can write Rd as a function of only two variables Rd = 
Rdicqz/cTxy, cr'^ /c). Thus, by similar arguments to what 
we had for the semi-infinite case, the Casimir pressure 
Pc = Pof{(7xya/c,al^/c). 

The limiting cases can be understood now by con¬ 
sidering first Eq. 0. The exponential constrains 
qz ~ 1 /a and since technically the thin-film limit is 
liInq^d-^o Rdiiqz) = 7 ?o(f( 7 z), we have that d/a —>■ 0. In 
other words, the thin film limit is applicable when we 
are considering d a. The opposite limit is just when 
qzd — 7 > 00 , and by similar arguments, that means d ^ a 
is when the semi-infinite case applies. Both limits leave 
a^ya/c and a^yd/c unaffected (although in the thin film 
case axya drops out whereas in the semi-infinite case 
axyd —>■ 00 has the same limit as qzd ^ 00 ). 

The thin film limit can be evaluated 

exactly [T5] and has the form = 

-Po|TRe{Li 4 [(cr^°/c) 2 /(cr 2 D/c-k 2 *) 2 ]} where LU is 
the polylogarithm of degree 4. Note that this function 
has a minimum value of ~ —0.117Po representing 
how repulsive we can get. For large enough cr™/c, the 
force does become attractive—corresponding roughly 
to when (a^y/c)^ > cr™/c (i.e. when Kerr rotation is 
suppressed). 

The cross-over between these regimes can be seen in 
Fig. ^ As Gxydjc is increased, the Casimir energy ap¬ 
proaches the semi-infinite case. However, for any finite d, 
each curve asymptotically approaches its thin-film value 
(and never goes lower than the minimum value repre¬ 
sented by the dashed horizontal line in Fig. |^. This not 



FIG. 3. (Color online) A plot of the normalized Casimir force 
for various thicknesses of a bulk Hall material (idealized Weyl 
semimetal). It begins slightly repulsive for small (Jxyd/c, and 
as this increases, it becomes more repulsive until it reaches the 
maximum for a thin him material (the dashed line) at which 
point it increases to the semi-inhnite limit. Pq = — and 

Oxy = e^b/2Ti^h is the bulk Hall response. Figure|^takes into 
account more material properties. 


only clearly connects our case to the previously known 
thin-film result, but also provides a theoretical justifi¬ 
cation for considering a thin-film limit d <C a with a 
two-dimensional conductance (T^,jd. 

Until now the plates have been idealized. Using the 
thin-film limit illustrated above as a reference allows us 
to easily consider some of the effects of taking into ac¬ 
count the full frequency response of the material. Thus, 
we pick a axyd that is reasonably in the repulsive regime 
(for all distances) in order to analyze the effects of includ¬ 
ing some of the lowest-order frequency dependence into 
the conductivities. We will mainly be interested in the 
effects of virtual vacuum transitions that are low in en¬ 
ergy, which correspond to plates that are far apart from 
one another. Thus, we will use the low-energy chiral 
Hamiltonian for a pair of Weyl nodes, 

Hy[r = HivFcr ■ (k ± b). 


( 8 ) 

















4 



FIG. 4. (Color online) The Casimir force for a thin film Weyl semimetal taking into account low-energy virtual transitions 

we 


in the band structure. An antitrap develops when the longitudinal conductivity overwhelms the Hall conductivity. In (a) 
compare different values of b (or equivalently, changing the Hall effect), in m we compare different up’s (the larger up, the 
smaller a^x is), and in (c) we turn on a finite chemical potential which causes attraction at very long distances (and hence 
a trap). Even small chemical potentials have this property, but the trap is quite far out. Unless the parameter is varying, 
ao = 1 nm, d = 20 nm, b = 0.3(27r/ao), A = 27r/oo, up = 6 x 10® m/s, and fj, — 0. 


where up is the Fermi velocity and b is the position of the 
of Weyl node in k space. The exact band structure will 
be important as the plates get closer although weighting 
will still be larger on the lower-energy modes. 

To the conductivities, we fix and calculate two- 
dimensional conductivities using the Kubo-Greenwood 
formulation (see the Supplemental Material [2H] for de¬ 
tails) then integrate the resulting expression over kz with 
a symmetric cutoff = j^^^a^^(iui\kz) [IS] 

where d^^{iuj]kz) is the two-dimensional conductivity 
with kz fixed. We evaluate this at an imaginary fre¬ 
quency to aid the Casimir calculations. 

We perform this procedure at hnite chemical potential 
/i and throw out terms that go to zero when the cutoff 
A —)• oo. This yields the conductivities. 


c(*w) = 


127r^fiz;F 


[|^ + 2^1og(^) 




log -I- , (9) 


and (Jxy{iuj) = 2 ^^ is unchanged at this order. Due 
to the linear dispersion of the Weyl nodes, we have a 
logarithmic cutoff dependence. Note that rotating to real 
frequencies we get the correct result for two Weyl nodes 
for Re[tTa;x(a;)] [SDj, and a result with the appropriate 
logarithmic divergence for ISSj- This can be 

understood in terms of charge renormalization due to the 
band structure, but for ease of our purposes we let A ^ 
1/ao where ag is the lattice spacing. For our plots we 
choose a lattice spacing of oq = 1 nm, a thickness of 
d = 20nm, b = 0.3(27r/ao), A = 2tt jag, wp = 6 x 10® m/s, 
and /r = 0 unless its the parameter we are varying. 

Now, one can use one of two equivalent ways of calcu¬ 
lating the Casimir energy: the reflection matrix as given 
in Ref. [T2|, or using a microscopic analysis to find the 
photon dressed RPA current-current correlators [28]. In 
order to avoid an unphysical negative as well as 

for consistency, we cutoff the photon energies in the Lif- 


shitz formula to run from 0 to A—an approximation valid 
for a » -^A“A 

Up 

First, we see that we get an anti-trap for these at ap¬ 
proximately 650 nm, and if we increase b as in Fig. [^ 
(with, say, an applied magnetic held), it not only moves 
closer to zero separation, but the overall repulsive behav¬ 
ior can be enhanced. On the other hand, if we increase 
Up as we see in Fig. [^ we see the region of attraction is 
suppressed, but the overall repulsive behavior at long dis¬ 
tances is maintained. Modifying A will have effects simi¬ 
lar to modifying up, but since it appears logarithmically, 
it would need to change by orders of magnitude to give 
appreciable changes (a simple plot for this is provided 
in the Supplemental Material [23 but is not relevant for 
the discussion here). This antitrap effect is occurring at 
short distances when higher order band effects also play 
a role, but any other effects will contribute to the longi¬ 
tudinal conductivity in such a way that an anti-trap will 
appear. 


Interestingly, when we introduce a finite chemical po¬ 
tential as we see in Fig. 4c in addition to the anti-trap 
we get at shorter distances, we start to see a trap at much 
longer distances appear. This is not surprising since at 
zero frequency there is a divergent longitudinal conduc¬ 
tivity. Thus, we know that at long distances, the Casimir 
force must be attractive, but by modifying the Hall ef¬ 
fect, we have an intermediate regime of repulsion. A 
similar effect would occur if we took finite temperatures 
or disorder corrections to the longitudinal conductivities. 


Considering the form of the conductance in terms of 
the fine-structure constant Uxydjc = we see that 

bd controls the strength of the repulsion in the thin-film 
limit. Without longitudinal conductance, the repulsive 
regime roughly corresponds to when {axydlc)"^ < Oxyd/c 
or equivalently ^ ^ The longitudinal conductance 
introduces vp into the scheme, relevant photons have w « 
c/a, and thus it becomes important for Oxxd/c 
0(1) (neglecting constants) which both emphasizes that 
Up controls the longitudinal conductance’s contribution 
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to the Casimir effect and that the term is suppressed at 
longer distances. 

We have shown here how Weyl semimetals can ex¬ 
hibit a tunable repulsive Casimir force (with, for in¬ 
stance, magnetic-field tuning b) and how it can depend 
on the thickness of the material. In the thin-film limit, 
we showed how the semimetallic nature of these mate¬ 
rials can work to create attraction at shorter distance 
scales and how a finite longitudinal conductivity will cre¬ 
ate long-distance attraction along with repulsion at in¬ 


termediate distances. Recently the first experimental ob¬ 
servation of Weyl semimetals jSl] [32] provided optimism 
that these theoretical materials could be a reality. The 
marginal nature of these materials could be useful for 
controlling the Casimir force between attractive and re¬ 
pulsive regimes. 

This work was supported by the DOE-BES (Grant 
No. DESC0001911) (A.A.A. and V.G.), the JQI-PFC 
(J.H.W.), and the Simons Foundation. We thank Liang 
Wu and Mehdi Kargarian for discussions. 
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SUPPLEMENTARY MATERIALS 
I. AXIONIC ELECTRODYNAMICS 

In the main text, we mention the axionic term which appears in the action alongside the usual Maxwell action 

Sa = J (10) 

where e is the electric charge, h is Planck’s constant, c is the speed of light, = d^Ai, — is the electromagnetic 
field tensor, is the fully antisymmetric 4-tensor, and 0(r, f) = 2b • r — 2bot is the axionic field. 
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For our purposes, we will set h = bz and bg = 0. If we only apply S'a for z > 0, there is no resulting surface current 
(i.e. this is the surface without Fermi arcs), and the current response is 


i’W = 
i“W = 

Now, we solve Maxwell’s equations in the bulk Hall system after taking the Fourier transform 

k • E = —iaxyZ ■ B, 
k • B = wB, 

k X E = wB, 
k X B = ia^yi x E — wE. 

One can define a frequency-dependent dielectric permitivity e{ui) to be 


e(cu) — I iOxyj^ 
0 


iOxyliO 0^ 

1 
0 


( 11 ) 

( 12 ) 


(13) 

(14) 

(15) 

(16) 


(17) 


in which case, Maxwell’s equations can be recast as a single equation for the electric displacement D = e(w)E, 

[k( 8 )k-fc 2 l]e-i(w)D = -a; 2 D. (18) 

The determinant of this matrix equation yields the frequencies that a wave vector k can have. In our case, we obtain 

(19) 


wit = -b ^^ly^(^xy\lk'i 


lcr2 

i xy 


There is a polarization associated with each of these frequencies which we can obtain from Eq. (18 1 using k • D = 0 
(no free charge). We choose ei = y x k and 62 = y as our basis for the polarizations (assuming ky = 0 without loss 
of generality). The resulting (unnormalized) polarizations are 


Di 2 — 


U}± 


kl 


1 2 1 
i'^xy T 2 '^xy 


ei ± ikze2. 


( 20 ) 


Notice that these polarizations are elliptical. 

To find the reflection of a wave off a half space hlled with this material, kx-, ky, and w must remain the same on 
either side of the material, but kz can change, and matching both sides of the dispersion Eq. (19), we obtain simply 
that for an incident wave with wave-vector q = iqx,Qy,Qz), the transmitted wave has 


{kz ) — Qz{qz i ^xy)- 


( 21 ) 


Each of these can be associated with a (unnormalized) polarization as (assuming qy = 0 without loss of generality) 

'D± = -^{qz±<Jxy)eiTikze2- ( 22 ) 

At this point, we note some interesting electromagnetic properties appearing here. Eq. ( |21[ ) implies that this material is 
birefringent, and for an incident wave (at any angle) with qz < \crxy\, only one (elliptical) polarization even propagates 
into the material while the other is an evanescent wave - independent of the angle of incidence. 


A. Reflection coefficient of semi-infinite bulk 

Now, to obtain the reflection matrix, we call our incident wave Eg with wave-vector q, our reflected wave E,. 
with wave-vector q^ = {qx,qy,—qz), as well as E± = e“^(w)D± with wave vectors k± for the two polarizations it is 
transmitted into. The relevant Maxwell equations at the interface between vacuum and the bulk Hall material are 
then given by 

(Eo-bE^-E+-E_) xz = 0, 

(q X Eg -f qr X Er — k 4 . X E_|. — k_ x E_) x z = 0. 
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We can break up the polarization of the incident and reflected waves into transverse electric (TE) and transverse 
magnetic (TM), and the reflection matrix is dehned such that 



Riu},q) 



Solving for this matrix, we obtain 


^oo(w,q) 


1 


^ xy 


( ^xy + k^ - k+ i{2q:, - - k+ )\ 

-k- -k+) a^y + k- - k+ ) 


This can then be rotated to imaginary frequencies and the result is in the main text. 


(24) 


(25) 


B. Reflection coefficient of thickness d sample 


If we have a material of thickness d, then we need additional matching conditions due to Maxwell’s equations at the 
other interface. This requires restricting our action Eq. ([^ to 0 < z < d. To solve this, we just need to add another 
set of matching conditions. In addition to the incident Eq and reflected waves, we now have forward moving Weyl 
polarizations E^ with = {kx,ky,kf), backwards moving Weyl polarizations E^ with = {kx, ky, —kf), and a 
transmitted wave Ej with wave-vector the same as the transmitted k. 

The resulting matching conditions are 




(Eo -b 

Er - 

-e;- 

eL - 

E:^ - 

Ei) 

X 

z = 

= 0, 

(k X Eo 

X 

1 

X 

-b 

— kl X 

e1 - 

- k);. X 

- 

ki X 

Ei) 

X 

z = 

= 0, 


(ETg*fc+d-^Ele*'=-‘* 

-bE)^e- 

-ifc d 

-bEie 

— ik~ d 

-Ete' 

ikzd'j 

X 

z = 

= 0, 

(k;x 

eT gifej-d kl X Ele*'”-'' 










-bk:^ X 

7<i^ki 

X E 

b g — *^7 


X Ete 

ikzd-^ 

X 

z = 

= 0. 


These equations can still be solved and the result is 


Rd{oJ,q) 


( Rxx 
\—Rxy 



(27) 


where 


Rxx = 0 ' 2 ,y{sin(fcJ(i)[jfc^ cos{k^d) + axySin{dk^)] — ik~ cos{k~d) sin{k^d)}/D, (28) 

Rxy = ^xy{kj cos{k~d) sin{k^d) + sin(A:“d)[fc^ cos(fc^ d) — 2ikz sm{k^d)]}/D, (29) 

D =[2ik~ cos(k~d) + (2kz — axy) siii(k~d)][2ik^ cos(fc+d) -|- (2kz -f axy) siii{k^d)]. (30) 


And again, this can be rotated to imaginary frequencies to obtain the result in the main text. 

And as stated in the text, the various limits (semi-infinite to thin film limits) apply to this reflection matrix. 


II. CALCULATING THE CONDUCTIVITIES 

In order to calculate the conductivities in the clean limit we consider the Hamiltonian near a Weyl node 

H\y = ihvFcr ■ (k ± b). (31) 

For simplicity we set h = 1 = vp unless otherwise specified. We consider a pair of these nodes, calculating the 
quantities separately for each node and adding them together (which will just introduce a factor of two for both axx 
and axy). 

To And the conductivities, we find a complete basis of states which are easily found by diagonalizing the two-by-two 
matrix Hw in momentum space (label them |/fc±)). The relevant matrix elements are then (choosing the negative 


sign for the Hamiltonian) 


{fk+Wx\fk-) 

{fk+\cry\fk-) 

{fk±Wx\fk±) 


{kz -h)kx+ ieky 

ex/e" - (fcz - 

{kz - h)ky - iekz; 

- (kz - ’ 



(32) 

(33) 

(34) 


where e = + ky + {kz — by. 

We then fix kz and use the Kubo-Greenwood formula for the intra- and inter-band transitions separately to obtain 
two-dimensional conductivities. Thus, 




6 ^k7 ^k7' (/k7|j^l/k7') (/k7'|/k7) 

^ k 7^7^ ^k7 ^k7' “t“ f^k7 ^k7' 


(35) 


where is the occupation in that band of that momentum and jy = are the single particle current operators. 
For intra-band quantities with nk± = 6*(/^ T Ck) 


-intra / ■ n \ '' ,, ^ X {fk+\jyfk+) (/k+|>|/k+) 

/ ^k,+ ) 5 

^ I “ ILU 


(36) 


assuming only the upper-band for simplicity and without loss of generality (due to particle-hole symmetry). 
Adding these contributions together at hnite chemical potential yields [T2] 


^XX (^^5 kz^ - 


a^yiiuj] kz) = 


o 2 r 


47r 

r 


1 - 


4:{kz -b)‘^\ i 


27r 


kz-bi . 

- o log 

uj 2 V 


log 


/ 2A — ioj 


^ 2 A — iuj 


2A 


A' 

ui 


2A 


- lUJ 


where A = max{|fcz — b\, |/r|}. With these quantities we can then use the cutoff procedure explained in 

I-A 


ay„{iU}) = / -^cry^(iw,kz). 

J-A 


(37) 

(38) 


(39) 


Through which we obtain (throwing away terms that go to zero as the cutoff increases to infinity and multiplying by 
two for the two nodes and bringing back in the constants h and fp) 


^XX (^^) 
(7xy (^^) 


Vl'K'^hvp 

e^b 

2Tr^h 


|a. + 2a;log(^)+4^-a;log(l + 


j\ ’ 


These quantites are what we use in the next section as input for the Casimir Force. 


(40) 

(41) 


III. CASIMIR FORCE CALCULATION 
A. Differentiating the Casimir energy 


In the main text, we find the Casimir pressure by taking the derivative of the energy = —dEc/da. This leads to 
the expression for force found in the text for two semi-infinite Weyl plates 


Pr = 


2hc 


dq, 


<ll 9 




2qza 


(42) 
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a (/im) 

FIG. 5. The Casimir energy for the conductivities as defined in the main text with b = 0.3(27r) nm~^, wf = 6 x 10® m/s, and 
fj, — 0. The cutoff is varied here. 


and the function g{u,v) is defined by 


, - [Rxxiu)'^ + Rly{u)]'^e 

- 2[RUur - Rlyiu)] - [RUu)^ + RlyiuWe-^ ’ 

where Rxx{u) and Rxy{u) are the matrix elements of Raoi^c^z ! xy) (Eq- (3) in the main text), 


Rxx(u) = \j2u{^/vF^- 1 ) - 1 , 
Rxy{u) = 2u — \j 2u{\/u‘^ + 1 + 1). 


(43) 


(44) 

(45) 


B. Two-dimensional plates calculation 

Another well-known approach completely equivalent to the Lifshitz formula comes directly from quantum field 
theory. 

To put it briefly, if we have a conductivity like we calculated above and write it as a response function n(za;) such 
that 


a{iuj) = —Il{iu!)/uj 


(46) 


And calculate the RPA response function considering photons “skimming” along the surface of our material 


n(iw, g) = [I — Il{iuj)D{iu!, q, 0)] ^ n(zw). 


(47) 


with the photon propagator 


D{tu;,q,z)={ 


(48) 




then the Casimir energy takes the form 


1 f 

Ec{a) = - / ^ / dwtrlog I - ni{iuj,q)D{iuj,q,a)Il 2 {iuj,q)D{iuj,q,a) 

TT J (27r)^ Jq 


(49) 


For our particular case of the conductivities introduced in Section [T^ we need to also cutoff the frequencies in this 
integral to go from 0 to upA for consistency. We expect higher energy virtual photons to not play a large role. 

We also show the result we obtain from varying the cutoff in Fig. Note that it does not affect the effect much 
unless it varies by orders of magnitude. 






















